Direction Dependence of Spin Relaxation in Confined 2D Systems 
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The dependence of spin relaxation on the direction of the quantum wire under Rashba and Dres- 
selhaus (linear and cubic) spin orbit coupling is studied. Comprising the dimensional reduction of 
the wire in the diffusive regime, the lowest spin relaxation and dephasing rates for (001) and (110) 
systems are found. The analysis of spin relaxation reduction is then extended to non-diffusive wires 
where it is shown that, in contrast to the theory of dimensional crossover from weak localization to 
weak antilocalization in diffusive wires, the relaxation due to cubic Dresselhaus spin orbit coupling 
is reduced and the linear part shifted with the number of transverse channels. 



PACS numbers: 72.10.Fk, 72.15.Rn, 73.20.Fz 



I. INTRODUCTION 

Spin dynamics in semiconductors have been studied for 
decades, but still the prime condition for building spin- 
tronic devices, namely the understanding of spin relax- 
ation, is not satisfactorily fulfilled. In the following we 
focus on materials where the dominant mechanism for 
spin relaxation is governed by the D'yakonov-Perel S pi n 
relaxation (DPR). 1 This mechanism results from lifting 
the spin degeneracy which is due to time inversion sym- 
metry and spacial inversion symmetry and leads to the 
effect of slower spin dephasing the faster the momen- 
tum relaxes (motional narrowing 2-4 ). In one of the most 
studied systems GaAs/AlGaAs DPR is the most relevant 
mechanism in the metallic regime. 5 

Preserving time reversion symmetry, the spin splitting 
can be due to bulk inversion asymmetry (BIA) 6 and also 
due to the asymmetry arising from the structure of the 
quantum well (QW), the structure inversion asymmetry 
(SIA) 7 . In Ref. 8,9 it was shown how the spin relaxes in 
a quasi ID electron system in a QW grown in the [001] 
direction, depending on the width of the wire, where the 
normal of the boundary was pointing in the [010] direc- 
tion. It is already known that in a (001) 2D system with 
BIA and SIA we get an anisotropic spin- relaxation. 10-12 
This has also been studied numerically in quasi- ID GaAs 
wires 13 . In this work, Sec. II, we present analytical results 
concerning this anisotropy for the 2D case as well as the 
case of QW with spin and charge conserving boundaries. 
We also extend our analysis to other growth direc- 
tions, Sec. Ill: Searching for long spin decoherence 
times at room temperature, the (110) QW attracted 
attention. 14,15 The properties of spin relaxation in sys- 
tems with this growth direction has also been related 
to weak localization (WL) measurements, Ref. 16. We 
present analytical explanations for dimensional spin re- 
laxation reduction and discuss the crossover from WL to 
weak antilocalization (WAL), Sec. IV. 



As we will show in the following sections, the cubic Dres- 
selhaus spin orbit coupling (SOC) gives always rise to 
a limitation of the spin relaxation time in the diffusive 
case, W S> l e , with the wire width W and the elastic 
mean free path l e . However some of the experiments are 
done on ballistic wires and we need to modify the theory 
used in Ref. 8,9 to enable us to study the crossover from 
diffusive to ballistic wires. In Sec. V we show how the 
spin relaxation which is due to cubic Dresselhaus SOC 
reduces with the number of channels in the QW. 
We consider the following Hamiltonian with SOC 
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H = ^t(p + eA ) + y W " (B + B so (p)) , (i) 



2m 



where 



is the effective electron mass. A is the vector 



potential due to the external magnetic field B. Bg = 
(B SOx , B SOy ) is the momentum dependent SO field, <x 
is a vector, with components <r%, i = x,y,z, the Pauli 
matrices, 7 is the gyromagnetic ratio with 7 = g/x B with 
the effective g factor of the material, and fx B = e/2m e is 
the Bohr magneton constant. For example, III-V and II- 
VI semiconductors such as GaAs, InSb have zinc-blend 
structure. This BIA causes a SO interaction, which, to 
lowest order in the wave vector k, is given by 6 
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where the principal crystal axes are given by i e 
{x, y, z}, i — > {{i — 1) mod 3) + 1 and the spin-orbit co- 
efficient for the bulk semiconductor . We consider the 
standard white-noise model for the impurity potential, 
V'(x), which vanishes on average (V(x)) = 0, is un- 
corrected, (V(x)V(x')) = <5(x — x')/27wr, and weak, 
cft S> 1. Here, v = m e /(2Tr) is the average density 
of states per spin channel, is the Fermi energy and 
t is the elastic scattering time. To address both, the 
WL corrections as well as the spin relaxation rates in the 
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system, we analyze the Cooperon 
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1 + zr(v(Q + 2eA + 2m e aS) + H a , + H z ) 



(3) 



where the integral is performed over all angles of velocity 
v on the Fermi surface, H a > = — (Q + 2eA)aer' and the 
Zeeman coupling to the external magnetic field yields, 



Hz = —jli**' - °") B - 



(4) 



The coupling between the orbital motion and the spin 
S = (er + <r')/2 is described by the SOC operator a. The 
spin quantum number is 1 instead of 1/2 due to the 
electron-hole excitation. It follows that for weak dis- 
order and without Zeeman coupling, the Cooperon de- 
pends only on the total momentum Q and the total 
spin S. Expanding the Cooperon to second order in 
(Q + 2eA + 2m e aS) and performing the angular inte- 
gral which is for 2D diffusion (elastic mean free path l e 
smaller than wire width W) continuous from to 2tt, 
yields: 



C(Q) 
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D e (Q + 2eA + 2eA s ) 2 + H. 



(5) 
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The effective vector potential due to SO interaction is 
As = m e aS/e, where a = (a) is averaged over angle. 
The SO term H JD , which cannot be rewritten as a vector 
potential, is in our case due to the appearance of cubic 
Dresselhaus SOC. 



with a 2 the Rashba parameter. 7,18 Therefore the 
Cooperon Hamiltonian, in the case of Rashba and lin. 
and cubic Dresselhaus SOC is given by 



H r := 



c- 



D, 



= (Q + 2eA s ) 2 + (m 2 e e FlD ) 2 (S 2 x + S 2 ), (8) 



with the effective vector potential 
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A * = V* 8 = V la? "a? OH S v I' W 



with 5i\ = otx — rrie'YD^F /2. 

It can be easily shown that the Hamiltonian Eq. (8) has 
only non vanishing eigenvalues due to (m 2 eFjD) 2 in the 
2D case. 

The term with (S 2 + S 2 ), which is due to cubic Dressel- 
haus SOC, is not reduced by reason of the boundary in 
the diffusive case. However two triplet eigenvalues of this 
term depend on the wire width, 
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(10) 
(11) 



with q 2 3 /2 = (m 2 eFlD) 2 ■ In the following we are go- 
ing to diagonalize the whole Hamiltonian and change the 
direction of the wire in the (001) plane. 



II. SPIN RELAXATION ANISOTROPY IN THE 
(001) SYSTEM 



A. Example 

To get an idea of the procedure we recall the situation 
presented in Ref. 8,9. Spin relaxation in a (001) quasi-lD 
wire in [100] direction: 

The Dresselhaus term, Eq. (2), for QWs grown in the 
[001] direction is given by 6 



-7B so . D 



ai(-e x p x + iyPy) +jD(e x PxPl ~ e y p y p 2 x 



(6) 

Here, ot\ = jjj(p 2 ) is the linear Dresselhaus parame- 
ter, which measures the strength of the term linear in 
momenta p x ,P v in the plane of the 2D electron system 
(2DES). When (p 2 z ) - I /a 2 > k 2 F (a is the thickness of 
the 2DES, kp, Fermi wave number), that term exceeds 
the cubic Dresselhaus terms which have coupling strength 
jd- Asymmetric confinement of the 2DES, a SIA, yields 
the Rashba term which does not depend on the growth 
direction 
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Ol2{e.xPy - e y Px), 



(7) 



A. 2D system 

We rotate the system in-plane through the angle 8 (the 
angle 9 = tt/4 is equivalent to [110]). This does not effect 
the Rashba term but changes the Dresselhaus one to 11,12 
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a v k v cos(20)«fc 2 ) - kl) - a x k x cos(20)((k 2 ) - k 2 ) 



-a y k x -sm(26)(k 2 x -kl-2(k 2 z )) 
+ a x k y ^sm(29)(k 2 x -k 2 y + 2(k 2 z )), 



(12) 



with the wave vectors The resulting Cooperon Hamil- 
tonian, including Rashba and Dresselhaus SOC, reads 
then 

H c = (Q x + a-xiSx + (a x2 - q2)S y ) 2 



+ (Q y + (a x2 + q 2 )S x - a xl S v Y 

+ q f(s 2 x + s 2 ), 



(13) 
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where we set 



9|a 
2 

CXxl 
&x2 



(m 2 e e F ~f D ) 2 



-m elD cos(29)((m e v) 2 - 4(k 2 z )), 
1 m elD sm(29){{m e v) 2 - A(k 2 z )) 



(14) 
(15) 
(16) 

(17) 
(18) 



with gi = 2m e ai, 52 = 2m e a.2- We see that the part of 
the Hamiltonian which cannot be written as a vector field 
and is due to cubic Dresselhaus SOC does not depend on 
the wire direction in the (001) plane. 




where n is the unit vector normal to the boundary dS 
and x is the coordinate along the wire. In order to do a 
diagonalization taking only the zero-mode into account, 
we have to simplify the boundary condition. A transfor- 
mation acting in the transverse direction is needed ac- 
cording to Eq. (13): C -)• C = U A CU\, by using the 
transformation 



U = 1 4 - ishi(q s y)—A y + (cos(q s y) - 1)\a. 

Is qj 



with A y = (a X 2 + q2)S x — a x \S y and 



q s = \J (a X 2 + <?2) 2 + ■ 



Spin relaxation 



(23) 



1. Special case: Only lin. Dresselhaus SOC equal to Rashba 

SOC 

As a special example for the 2D case we set q S 3 = 
and q\ =92- To simplify the search for vanishing spin 
relaxation we go to polar coordinates. Applying free wave 
functions (with k x ,k y ) to H c , Eq. (13), we end up with 
(singlet part left out) 



He 
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c.c. 




(19) 



with k x /q2 = Q cos ((/>), kyjq-x = Qsin((/>) and 
j H = (i-l)V2exp(id)(cos(<P + 6)-sm(<P + 9))Q. (20) 

Vanishing spin relaxation is found at Q = for arbitrary 
values of 9 (the spin with vanishing spin relaxation is 
pointing along the [110] direction 19 ). Another solution is 
found at Q = 2 with the condition 9 + (j> = 37r/4, which 
is equivalent to the [110] crystallographic direction. 11 

B. Quasi- ID wire 

In the following we consider spin and charge conserv- 
ing boundaries. Due to the SOC we have the following 
modified Neumann condition 9 



-— n • (v F [ 7 B so (k) • S]) - td n ) C\ ds = 0, (21) 



where (...) denotes the average over the direction of vj? 
and k which we rewrite for the rotated x-y system 

(-id y + 2e(A s UC (x 7 y = ±%- ) = 0, Wx, 



(22) 



We diagonalize the Hamiltonian, Eq. (13), after apply- 
ing the transformation U, taking only the lowest mode 
into account. The spectrum of the Hamiltonian for small 
wire width, Wq s < 1, is given by 



E 1/2 (k x > 0) 
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w 2 , 



(24) 
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E 1 (k x = 0)=qi 3 +q i 

v2 _i_ „,2 „2\2 , gs3 2 „2 



Ki + ^2-g|)^ + ^g; ^ 2 
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E 2 (k x = Q) = ^f + q 2 rn + W 2 



2 3q 
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Sill 



o q 2 (2 + Ism) { a ll + a l2 - qf> 



sm 



(25) 
(26) 

1 

-W 2 , 
(27) 



with q sm = y/ (a X 2 — q2) 2 + o? x \- First we notice that the 
only 9 dependence is in the term q sm , which disappears 
if the Dresselhaus SOC strength Si, which is shifted due 
to the cubic term, equals the Rashba SOC strength 0:2 
and the angle of the boundary is 9 = (1/4 + n)ir, n € 7L. 
Assuming the term proportional to W 2 /k x to be small, 
the absolute minimum can be found at 
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(28) 



4 



which is independent of the width W if a x \{9 = 0) = —q 2 
and/or the direction of the wire is pointing in 



e = 1 arcsin f 2{kl){ melD f{{m e vf-2{kl))-ql 
2 \ (mfu 2 7D - ^{k 2 )m e ^ D )q 2 

. (29) 

The second possible absolute minimum, which dominates 
for sufficient small width W and q sm ^ (compare with 
E2{k x = 0)), is found at 



E, ■ 
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q 2 

Is in 



l xl 
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x2 
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2 

sin 
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(30) 



2. Spin dephasing 

Concerning spintronic devices it is interesting to know 
how an ensemble of spins initially oriented along the [001] 
direction dephases in the wire of different orientation 6. 
To do this analysis we only have to know that the eigen- 
vector for the eigenvalue E\ at k. x = 0, Eq. (25), is the 
triplet state \S = l;m = 0) = (| U) + I it))/\/2 = 
| =£) = (0, l,0) r . This is equal to the z-component of 
the spin density whose evolution is described by the spin 
diffusion equation. 9 As an example we assume the case 
where cubic Dresselhaus term can be neglected and where 
the Rashba and lin. Dresselhaus SOC are equal. We no- 
tice that the dephasing is than width independent. At 
definite angles the dephasing time diverges - as for the 
in-plane polarized states with eigenvalue E 2 (k x = 0) - , 



The minimal spin-relaxation rate is found by analyzing 
the prefactor of W 2 in Eq. (30), Fig. (1). We see imme- 
diately that in the case of vanishing cubic Dresselhaus or 
in the case where a x \{9 = 0) = — q 2 we have no direc- 
tion dependence of the minimal spin relaxation. Notice 
the shift of the absolute minimum away from q\ = q 2 
due to q S 3 ^ 0. In the case of qi < (q s ^/\/2) we 
find the minimum at 8 = (1/4 + n)ir, n G TL, else at 
6 = (3/4 + n)ir, n € *Z, which is indicated by the dashed 
line in Fig. (1). 
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Figure 1: (Color online) Dependence of the W 2 coefficient in 
Eq. (30) on the lateral rotation {&). The absolute minimum 
is found for a x i(0 = 0) = —52 (here: qi/q2 = 1.63) and for 
different SO strength we find the minimum at 9 = (1/4 + 
% if qi < {qs3/V2) (dashed line: q x = (qs3/y/2)) 
and at = (3/4 + n)n, n £ TL else. Here we set q 3 3 = 0.9. 
The scaling is arbitrary. 



2D e q 2 {l - sin(20)) 



(31) 



which is plotted in Fig. (2). We have longest spin de- 
phasing time at = (1/4 + n)7r, n € TL. For 6 = 
(3/4 + n)7T, n G TL we get the 2D result T 2 = l/(4gfD e ), 
which is given by the eigenvalue of the spin relaxation 
tensor 9,20 ' 21 , 
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r 7 2 ((B so (k) 2 )% - (B so (k) iJ B so (k) j )) (32) 



to the triplet state \S = 1; m = 0). 

This gives an analytical description of numerical calcula- 
tion done by J.Liu et at, Ref. 13. 

Switching on cubic Dresselhaus SOC leads to finite spin 
dephasing time for all angles 8. In addition T 2 is than 
width dependent. In the case of strong cubic Dresselhaus 
SOC where g 2 3 /2 = qf = q 2 , the dephasing time T 2 is an- 
gle independent and for q 2 3 /2 > qf = q\ the minima in 
T 2 (0) change to maxima and vice versa. 



3. Special case: 8 — 

In this case the longitudinal direction of the wire is 
[100]. 

If we neglect the term proportional to W 2 /k x in Eq. (24) 
the lowest spin relaxation is found to be 
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depending whether 
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Figure 2: (Color online) The spin dephasing time T 2 of a spin 
initially oriented along the [001] direction in units of (-De?!) 
for the special case of equal Rashba and lin. Dresselhaus 
SOC. The different curves show different strength of cubic 
Dresselhaus in units of q s i/q 2 . In the case of finite cubic 
Dresselhaus SOC we set W = 0A/q 2 . If q s3 = 0: T 2 di- 
verges at = (1/4 + n)n, n € 7L (dashed vertical lines). The 
horizontal dashed line indicated the 2D spin dephasing time, 
T 2 = l/(4g 2 2 D e ). 



we end up with the following Cooperon Hamiltonian 

= (Qx - QlS z - q 2 S v f + (Q y +q 2 S x ) 2 + ^S 2 z . (39) 

with qx = 2m e f (fc2) - ^^f^, q 2 = 2m e a 2 and q 3 = 
(3m e ej?(jD /%))■ We see immediately that in the 2D case 
states polarized in the z-direction have vanishing spin 
relaxation as long as we have no Rashba SOC. Compared 
with the (001) system the constant term due to cubic 
Dresselhaus does not mix spin directions. Here we set 
the appropriate Neumann boundary condition as follows: 

(-idy + 2m e a 2 S x )C (x, y = ±yj = 0, Wx. (40) 

The presence of Rashba SOC adds a vector potential pro- 
portional to S x . Applying a non-abelian gauge transfor- 
mation as before to simplify the boundary condition, we 
diagonalize the transformed Hamiltonian (App. (Al)) up 
to second order in q 2 W in the 0-mode approximation. 



is negative or positive. This shows that the cubic Dres- 
selhaus term adds not only to the relaxation rate by a 
constant term but is also width dependent. However, 
this width dependence does not reduce the spin relax- 
ation rate below q 2 3 /2. 



III. SPIN RELAXATION IN QUASI-1D WIRE 
WITH [110] GROWTH DIRECTION 

To get the spin- relaxation in a [110] QW with Rashba 
and Dresselhaus SOC again we have to rotate the spacial 
coordinate system of the Dresselhaus Hamiltonian Eq.(2) 
but now with the rotation matrix 



A. Special case: without cubic Dresselhaus SOC 



J- ^~ 
1 



(36) 



We get 



H 



D[no] 
Id 



(t x [ k x k z 2kyk z + fc z ) 
+ a y (Ak x k y k z ) 

+ a z {kl~2k x k 2 y -k x k 2 z ). (37) 

The confinement in z-direction (z =[110]) leads to (k z ) = 
{kfj = 0, and (kf) = J \V(j}\ 2 dz. The Hamiltonian for the 
QW in [110] direction has then the following form 16 

ff[iio] = -lDo z k x Q<*2) - \{k 2 x - 2k 2 y )^j . (38) 

Including the Rashba SOC (q 2 ), noting that its Hamil- 
tonian does not depend on the orientation of the wire, 16 



The spectrum is found to be 

Ex = k 2 x + ±A 2 (q 2 W) 2 , 

E 2 , 3 = k 2 x + ±A 2 (2A-(q 2 W) 2 ) 
A 



(41) 



± —y/AifaW)* + 4fc?(24 - (q 2 W) 2 ) 2 , (42) 

with the lowest spin relaxation rate found at finite wave 
vectors fc, min = ±4(24 - (q 2 W) 2 ), 



1 A 2 



(43) 



We set A = ^q 2 + q\. 



B. With cubic Dresselhaus SOC 

If cubic Dresselhaus SOC cannot be neglected, the 
absolute minimum of spin relaxation can also shift to 
k Xmin = 0. This depends on the ratio of Rashba and lin. 
Dresselhaus SOC: 

If q 2 /qi <C 1, we find the absolute minimum at k Xmla = 0, 

E mml = g ~ 3 + f +ql A c + ±A c (q 2 W) 2 , (44) 
with 



= 2 \/(<?3 + <??) 2 + 2 (9i - I3)l 2 + 4- (45) 



(i 



If 52/91 ^> 1 , we find the absolute minimum at 
k Xnlin *±±(24~(q 2 W) 2 ), 



g2 

E m i n 2 = fc Xmin — fcx min <Z2 ( — + 2 
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/Ml 93 92 92 /- ~2 



V 12 3072fc3 mm 24 



(qs-Qt) 
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93 , 939? A 93 9| 



128 192 J 96 



W 2 . (46) 



We can conclude that reducing wire width W will not 
cancel the contribution due to cubic Dresselhaus SOC to 
the spin relaxation rate. 



IV. WEAK LOCALIZATION 

In Ref. 8,9 the crossover from WL to WAL due to 
change of wire width and SOC strength was explained 
in the case of a (001) system. Whether WL or WAL is 
present depends on the suppression of the triplet modes 
of the Cooperon. The suppression in turn is domi- 
nated by the absolute minimum of the spectrum of the 
Cooperon Hamiltonian H c . The findings presented in 
Sec. II B therefore point out that e.g. the crossover width, 
at which the system changes from WL to WAL, can shift 
with the wire direction 9. Recently experimental results 
on WL/WAL by J. Nitta et ai, Ref. 22, seem to show a 
strong dependence on growth direction. 
In the (110) system the situation is different: In the 2D 
case it was shown by Pikus et ai, Ref. 16, that in the ab- 
sence of the Rashba terms the negative magnetoconduc- 
tivity cannot be observed. In the case of a wire geometry 
we can conclude from Eqs. (41-46) that we have no width 
dependence if Rashba SOC vanishes. A change of the 
quantum correction to the static conductivity therefore 
cannot be achieved in this wire geometry by changing 
the wire width. The reason is the vector potential in the 
boundary condition, Eq. (40), which only depends on the 
Rashba SOC. 



V. DIFFUSIVE-BALLISTIC CROSSOVER 

In the following we assume a (001) 2D system with 
both, Rashba and linear and cubic Dresselhaus SO cou- 
pling. 

Experiments measuring WL in diffusive QW with 
SOC 23,24 are in great agreement with theoretical calcula- 
tions by S. Kettemann, Ref. 8. But considering e.g. the 
works Ref. 25,26, one realizes that the scope of applica- 
tion of the theory has to be extended to describe also the 



crossover to the ballistic regime, l e > W. We have shown 
in Sec. II B that the presence of cubic Dresselhaus SOC in 
the sample leads to a finite spin relaxation even for wire 
widths Q so W <C 1, regardless of the boundary direction 
in a (001) system. To account for the ballistic case we 
have to modify the derivation of the Cooperon Hamilto- 
nian, Eq. (8). In the case of a wire where the mean free 
path l e is comparable to the wire width W, we cannot in- 
tegrate in Eq. (3) over the Fermi surface in a continuous 
way. Instead, we assume kp/W to be finite and sum over 
the number of discrete channels N = [kpW/ir], where 
[. . .] is the integer part. Because H 1 d ~ e 2 F this constant 
term due to cubic Dresselhaus should reduce if we reduce 
the number of channels. If we expand the Cooperon to 
second order in (Q + 2eA + 2?7i, e aS) before averaging over 
the Fermi surface, (...), and use the Matsubara trick, we 
get 



^- = 2/i (o y + 2a 2 S x + 2 ( ai - 1dv 2 j^ 1 



.2/3 



2/2 ( Qx - 2a 2 S y - 2 I ai - 7d^J S 



+ 8j 2 D v 4 



(47) 



with m e = 1 and functions fi(tp) (App. B) which depend 
on the number of transverse modes N. In the diffusive 
case we can perform the continuous sum over the angle 
tp in Eq. (B3)-(B7), and we receive the old result with 
h = h = 1/2, h = 1/8 and h = h = 1/16: 



H c = {Q y + 2a 2 S x + 2 ( ai - -j D e F ) S v f 



1 



+ (Q x - 2a 2 S y - 2 yan - -jd^f ) S x ) 
+ ( lD e F ) 2 {S 2 x + Sl). 



A. Spin Relaxation at QsoW <JC 1 



(48) 



In the first section we analyzed the lowest spin relax- 
ation in wires of different direction in a (001) system. We 
have shown, that for every direction there is still a finite 
spin relaxation at wire width which fulfill the condition 
Q so W <C 1 due to cubic Dresselhaus SOC. It is clear 
that this finite spin relaxation vanishes when the width 
is equal to the Fermi wave length A^. In the following 
we show how this finite spin relaxation depends on the 
number of transverse channels N . We show in Ref. 27 
that the findings are consistent with calculations going 
beyond the perturbative ansatz. This is possible in a 
similar manner as has been done previously in Ref. 28 for 
wires without SOC, which showed the crossover of the 
magnetic phase shifting rate, which had been known be- 
fore in the diffusive and ballistic limit, only. 
To find the spectrum of the Cooperon Hamiltonian with 



7 



boundary conditions as in Sec. (II B), we stay in the 0- 
mode approximation in the Q space and proceed as be- 
fore: According to Eq. (47), the non-Abelian gauge trans- 
formation for the transversal direction y is given by 



U = cxp 



2a 2 S x + 2 I ai - "Idv 2 ^- ) 



fi 



(49) 

To concentrate on the constant width independent part 
of the spectrum we extract the absolute minimum at 
Q = 0, Fig. (3) and Fig. (4). A clear reduction of the 
absolute minimum is visible. Due to the factor /3//1 in 
the transformation U, the decrease of the minimal spin 
relaxation depends also on the ratio of Rashba and linear 
Dresselhaus SOC. 




Figure 3: (Color online) The lowest eigenvalues of the confined 
Cooperon Hamiltonian Eq. (47), equivalent to the lowest spin 
relaxation rate, are shown for Q — for different number of 
modes N = IcfW/tt. Different curves correspond to different 
values of a.2/q B - 




Figure 4: (Color online) The lowest eigenvalues of the confined 
Cooperon Hamiltonian Eq. (47), equivalent to the lowest spin 
relaxation rate, are shown for Q — for different number of 
modes N = IcfW/tt. Different curves correspond to different 
values of ai/q„. 



From Eq. (47) it is clear, that not only the H lD is af- 
fected by the reduction of the number of channels N but 



also the shift of the lin. Dresselhaus SOC, ot\. in the or- 
bital part. A model to extract the ratio of Rashba and 
lin. Dresselhaus SOC developed in Ref. 29 by Scheid et 
al. did not show much difference between the strict ID 
case and the non-diffusive case with wire of finite width. 
The results presented here should allow for extending the 
model to finite cubic Dresselhaus SOC. Deducing from 
our theory, the direction of the SO field should change 
with the number of channels due to the mentioned N 
dependent shift. 



VI. CONCLUSIONS 



Summarizing the results, we have characterized the 
anisotropy and width dependence of spin relaxation in 
a (001) QW. There are special angles 9 which are opti- 
mal for spin transport in quantum wires of finite width: 
The [110] and the [TlO] direction. At [110] we find the 
the longest spin dephasing time T 2 . If the absolute min- 
imum of spin relaxation is found at [110] or [110] direc- 
tion depends on the strength of cubic Dresselhaus and 
wire width. The findings for the spin dephasing time 
are in agreement with numerical results. The analytical 
expression for T 2 allows to see directly the interplay be- 
tween the cubic Dresselhaus SOC and the dimensional 
reduction, having effect on T 2 . In addition we analyzed 
the special case of a (110) system and found the mini- 
mal spin relaxation rates depending on Rashba and lin. 
and cubic Dresselhaus SOC in the presence of bound- 
aries. This results can be used to understand width and 
direction dependent WL measurements in QWs. Finally, 
we have shown how the reduction of channels in the wire 
reduces the finite spin relaxation rate which is due to cu- 
bic Dresselhaus SOC and does not reduce if the wire is 
small, Wq s <C 1, and diffusive, W ^> l e . The change in 
channel number also changes the shift of lin. Dresselhaus 
SOC strength, a\. This has to be considered if extract- 
ing SOC strength from wires with only few transverse 
channels. 
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Appendix A: Hamiltonian in [110] growth direction with wave vector Q. We set 



The Cooperon Hamiltonian in the 0-mode approxima- 
tion is given as follows 



ABC 
B* D E 

C* E* F 



M, 



q3 ■ 



(Al) 



with 

A = 



Aq 2 W 



(q 2 (Ak 2 x + 3(q 2 +q 2 ))W 



16fc x gi sin {^Y^j + (<7i - ql) sin(q 2 W)), 



B = 
C = 
D = 

E = 
F = 



lUk x sm(^)-q lS m(q 2 W) 



V2W 

_ Q2 (gf + gf) W + (q 2 - qj) sin(q 2 W) 
4q 2 W 

q 2 (2fc 2 + q 2 + qj) W + (q 2 - qf) sm{q 2 W) 



2q 2 W 

(4k x sm(&^}+q lS m(q 2 W) 

((ft (Ak 2 x + 3(q 2 +q 2 ))W 
q 2 W 



(A2) 
(A3) 
(A4) 
(A5) 

(A6) 



1 



Aq 2 W 

lQk x q\ sin 



q( - q$) sm(q 2 W)) (A7) 



and the term due to cubic Dresselhaus SOC 



M, 



q3 



\smc(q 2 W) + l \smc{q 2 W)-\ 
q 3 | \-\ sinc(g 2 W0 



isinc(g 2 W)-i 



\smc(q 2 W) 



(A8) 



Appendix B: Summation over the Fermi Surface 

The Cooperon Hamiltonian in the 2D case is given by 

H c = tv 2 {(cos 2 (ip)) (Q + 2m e a.S) 2 x 
+ (sin 2 (^))(Q + 2m e a.S) 2 
+ 4m 2 e j D v 2 (cos 2 (<p) sin 2 (v3))(Q + 2m e a.S) x .S x 
- 4m 2 7DW 2 (sin 2 (^) cos 2 ((^))(Q + 2m e a.S) y .S y 
+ (2ml lD v 2 ) 2 ((cos 2 { V ) sm 4 ^))S 2 



+ (sm 2 (p)cos i (p))S 2 )}, 



(Bl) 





= h 


(B2) 


fi 


= ( sin2 M), 


(B3) 


h 


= (cos 2 ((^)), 


(B4) 


h 


= (sin 2 (^)cos 2 (^)), 


(B5) 


fi 


= (sin (if) cos 2 (tp)), 


(B6) 


h 


= (sm 2 (ip) cos 4 (<p)). 


(B7) 



Using the Matsubara trick we write 



27r dip _ 2 ^ 

2tt ~~ TTiV ^ 



v 1 ^ 



This gives us 



h = 

h = 
h = 



2 N ~ X 



, N j 



N 



n N E \/ 1 ( N 



N 



2 ±(±)\, 



ttN ^ \N 

s=l 
N 



ttN ^ \N 

s=l 



. 3 
2 \ 2 



(B8) 



(B9) 

(BIO) 
(Bll) 
(B12) 
(B13) 



Writing Eq. (Bl) in a compact way gives us Eq. (47). 
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